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ABSTRACT 


A sea swell and surf model is implemented, tested and 
evaluated on a micro computer (HP-9845B). The two-dimensional 
model includes spectral wind-wave generation in open water, 
and shallow water wave transformation over irregular 
topography. The model predicts surf zone width, breaker 
lines and types of breakers. Using change in momentum flux 
of gravity waves (radiation stress) as forcing, the model 
predicts current velocities within the surf zone. The model 
is evaluated for the conditions over constant depths and 
uniform sloping beaches. The numerical results are checked 
against accepted theory and field observations. The model is 
found to overbuild wind generated sea heights for a 30 kts 
wind but to give expected wave heights for a 15 kts wind. 

The model results compare well with observed nearshore wave 
heights but give poor location of breaking waves. The 
model's nearshore current calculation is found unsuitable for 
the HP-9845B due to computational instability and time 
requirements. The sources of the model's problems are 
identified, and recommendations are made for future improve- 


ments of the model. 
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I. INTRODUCTION 


A. MOTIVATION 

Wave and surf forecasting 1S a primary concern of the 
U.S. Navy in executing amphibious operations. The wave 
environment affects the larger vessels in deep water and 
landing craft approaching the beach. As waves propagate from 
deep to shallow water, they undergo shoali ;, refraction, and 
eventually breaking, while forming a surf zone. The height 
and type of breakers in the surf zone determines the difficulty 
that landing craft will have in reaching the beach. The 
breaking waves can also generate nearshore currents; these 
shear currents potentially can turn a craft and cause it to 
be broached and swamped by breaking waves. In addition to 
amphibious units, coastal structures and pierside ships are 
subject to forces by waves in shallow water regions. 

The forecasting of wave conditions can be divided into 
two regimes: the nearshore, including the surf zone, and an 
offshore area. There are a number of numerical models avail- 
able which predict waves, shear currents and surf conditions 
separately. Propagating wave energy from deep water to the 
high water line on the beach requires a combination of comple- 
mentary models. The objective of this thesis is to convert a 
state-of-the-art computer program for use on a Hewlett-Packard 


9845B mini computer, to test the model against field observations 


and to revise where necessary. The program will be submitted 
as a prototype wave and surf model for shipboard use in the 


Fleet. 


B. BACKGROUND OF WAVE AND SURF PREDICTION 

The earliest models of wave generation and shallow water 
transformations were for monochromatic regular waves. Refrac- 
tion was based on the wave ray method (Arthur, 1949), where 
ray paths were calculated for individual waves with slowly 
varying wave speeds along the wave front. This method is 
analogous to ray path prediction used in optics for light 
rays. With the development of computers and a better under- 
standing of the physical processes, it became possible to 
consider the entire wave spectrum in generation and shallow 
water transformation. Karlsson (1969) introduced a direction 
variable to study wave spectral transformation for straight 
and parallel contours. Collins (1972) included bottom frveuied 
and evaluated the validity of one dimensional models. Noda, 
et al. (1974) used a relaxation finite-difference scheme to 
solve the stationary wave spectral transformation in shallow 
water with the inclusion of bottom dissipation and effects of 
local wind generation. Shiau and Wang (1977) studied wave 
transformation in shallow water and compared the results with 
field observations. Wang and Yang (1981) applied their model 
with the inclusion of bottom friction. Recently, Chen and 
Wang (1983) improved the model by considering non-stationary 


waves and interactions with currents. 


C. CURRENT STATUS OF MODELING TECHNIQUE 

The current state-of-the-art for numerical modeling can 
be characterized by two numerical methods, finite difference 
and finite elements. Both methods have been applied to wave 
and surf problems. A finite element model, such as by Wu and 
Liu (1985), economically handles localized irregular topography 
and has good properties for nonlinear convergence. However, 
finite element schemes are highly complex and require a sophis- 
ticated understanding of the model chara-teristics to choose 
the best grid system. Finite difference models (Noda, 1974; 
Shiau and Wang, 1977) use rectangular grids in general, and 
the numerical analysis has been widely studied. The numerical 
procedure of finite difference models is self-explained. There- 
fore, a finite difference scheme is less demanding on the 
user's numerical skills and from an operational point of view, 


is more practical. 


D. MODEL OVERVIEW 

The model described below was developed by Wang and Chen 
(1983). It is a finite difference model which considers a 
wind generated spectrum, spectral and swell wave transforma- 
tions over moderately irregular bottom topography, and 
dissipation by breaking waves and bottom friction. The model 
calculates a wind generated and deep water wave spectrum in 
the ocean and predicts wave and surf conditions on the beach. 


There is also an option of direct input of a predetermined 


wave into a surf area. In the nearshore computations, th:re 
are options for computing surf conditions and/or nearshore 
Circulation. The surf and nearshore circulation output are 
the basis for predicting amphibious landing conditions. Model 
input includes wind direction and speed, initial ocean swell 
and dominating direction, bottom type, and water depths. All 
the inputs are entered for the deep water calculation and 


again serve aS appropriate data for the shallow water. 


E.~ OBJECTIVES 

The objective of this thesis is to compile the numerical 
model by Wang and Chen (1983) in a user friendly version and 
test it on the HP-9845B mini computer. The program should be 
self-explanatory with internal error checks and simple correc- 
tion steps. The model will be tested against observations on 
several beaches to establish accuracy and range of application. 
The physical and numerical limitations of the model will be 
documented, explained and where possible improved. It should 


be a useful tool for the Shipboard Aerographer Mates. 
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Pie THEORY AND STRUCTURE OF THE MODEL 


A. THE SEA SWELL AND SURF PROGRAM (SSSP) 

The sea, swell and surf program (SSSP) was developed by 
Wang and Chen (1983) and is a combination of several comple- 
mentary models. The sea and swell for open sea follows the 
models by Shiau and Wang (1977) and Chen and Wang (1982). The 
nearshore current calculation is basec on a model by Birkemeier 
and Dalrymple (1976). 

The version of the SSSP most suitable for use on the HP- 
9845B consists of a main program which controls nine subrou- 
tines. A flow chart of the program is given in Figure l. The 
subroutine format of the SSSP allows for easier test and 
evaluation by component, modification or program correction, 
and reduced memory requirements on the mini computer. The 
main disadvantage of the subroutine format is some reduction 


in calculation speed. 


B. REFRACTION OF SEA SWELL 
The wave direction is determined using conservation of 


wavenumber for a steady state process (Phillips, 1977): 


tke cogue)  ouk sino) 


Jy Ox cane (1) 


where (x,y) are the on-offshore and alongshore directions, 


and 86 is the counter clockwise angle from the X axis to the 
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direction of swell propagation and k is the wavenumber. It 
is assumed that (1) is valid in both open water and within 
the surf zone. 

Noda (1972) developed an iterative scheme to solve (1) 
for the wave angle. The convention for grid location of the 
numerical model is shown in Figure 2. The differential equa- 
tion is central-differenced in the X-direction and forwarded- 
differenced in the Y-direction. The finite difference form of 
the equation usec by Noda (1972) with weighting on the forward 
difference term added to increase stability (after Abbot, 


1979) is used here: 


_ ee lire 
oF = sin im) $i1n O) i41,j-1 
pa 
(2) 
- AX ((, COS eo) - (k cos 86) yl 
2Ay ial ie, =a 
where weighting factor of t = 0.25 is used. The initial guess 


for the wave angle is calculated by applying Snell's law for 


parallel bottom contours. (This gives an exact solution in 

the case of a plaresloping beach.) The initial guess field 

Oe OS : is iterated until 05 : is within a half percent of 
rd f 


the previous iteration. 
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Figure 2. Differencing scheme used in the model. The 
grid points are identified by subscripts with 
tne first subscript indicating position from 
the beach with increasing numbers farther 
offshore. The second subscript indicates 
position alongshore with increasing numbers to 
the right when the beach is viewed from a 
ship off the coast. 
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The wave number ks : is the positive root of the linear 


a 


wave dispersion relationship 
¢ = (gktanh ka)l/? (3) 


where d is the water depth and o is the angular frequency. 

In the original SSSP model by Wang and Chen (1983), the wave- 
number was solved for in (3) using an iterative Newton's method. 
However, to speed up the model, the wavenumber is calculated 
here using a 6th degree polynomial fit of (3) as given by 
font (1979): 


p=" + (10 0.666675 + oms5550S- 


ga 


3 


+ 0.16084S° + 0.063208" + 0.02174§> 


(4) 


+ 0.006548" + 0.00171s’ + 0.00039s° 


+ 0.000118") ~1) 71 


where: 


57 (a/q) 


~” 
II 


This non-iterative approximation resulted in an observed 90 
percent savings in computation time on the HP-9845B with the 


same accuracy. 
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C. WAVE HEIGHT TRANSFORMATION 

The information of wave direction is used to calculate 
an initial guess for the energy field. The initial energyeae 
each grid point is calculated by assuming conservation of 


energy over straight and parallel contours. 


d0(E C_ cos 6) 
ern Ss 


ae = (5) 


Integrating (5) leads to the initial energy field: 


C cos @ 
O 


o C_ cos @ eal 
where the subscript (o) refers to conditions at the offshore 
boundary and 2 is the group velocity. The initial guess 
energy field is used to calculate the wave height field. 

The calculation of wave height is based on the energy 
balance equation including refraction, dissipation and growth. 
Three assumptions are made to speed calculation and reduce 
model complexity (Chen and Wang, 1983). First, it 1s assumed 
that wave energy within a particular spectral band stays in 
that band, which allows linear superposition of wavelets. The 
second assumption is that each frequency component can be 
described by a single mean direction, 6. Third, it is assumed 
that there are no wave-current interaction. With these 
assumptions, the steady state energy flux equation for each 
spectral band can be written (Longuet-Higgins, and Stewart, 
ng60 , 1962). 
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pg S—(S(£) C, cos 8) +o g ay (St£IC, sin 6) 


= - Eg (f) + €,(f) (7) 


where: 


S(f) = wave energy density at spectral frequency f 
o = density of the fluid 
Gees umecceteration due tO gravity 
€q(f) = energy dissipation through bottom friction 
E, (f) = energy generation due to wind stress 


The energy generation term, oat is based on the mechanisms 
used for wind energy transferred to the wave field described 
Pee nitlips (1957) and Miles (1957;1959a,b;1962). The Phillips 
mechanism is due to random atmospheric pressure fluctuations 
advected by wind resonantly interacting with a perturbed sea 
surface. The wave energy generated by the Phillips wind 


mechanism is 


S = QL (8) 

D OP g 
where 

An? ae 

a = 73 P(k,o) (3) 
p g 

and 

a, = a/g (10) 


fe 


Barnett (1968) fitted the following empirical relationship 


for P based on field measurements: 


_ 6.13 = Omens Me 
ERS) > > ee z 
TO v5 + (K Sino) 
v 
vi + (kK cos 67 — fy 
where: 
6 = angle between wind and wave 
w = wind speed (meters/second) 
A = o/w 
vy = 0.3 28 
7 035 
vee) = O52 


This mechanism is important for wave initiation from calm 
water and results ina linear increase in energy with time 
(Hasselman, 1960). 

The Miles energy transfer mechanism is important only 
after waves have formed (LaBlond and Mysak, 1978). When wind 
blows over an irregular sea surface, flow separation can 
occur in the lee of the wave crest creating pressure differ- 
ences. Momentum is transferred to the waves as a result of 
these pressure forces. The instability mechanism of Miles' 
results in exponential wave growth (Barnett (1968); Hasselmann 


C60) )o 
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S = og8(S(f) ] (12) 


m 
where: 
_ Aso,W cos 6 _ 
8 = Se z 5] CEs) 
s = ratio of air to water density 
A = 5 
B = 2 


The wind wave generation is then the sum of Sa) and (S_): 


€ = § +5 (14) 


The wave growth is limited to a fully risen sea at each 
spectral band. The limits on development are based on Phillips 
meo7) equilibrium spectrum modified by Kitaigorodski1 et al 


(1975) and Thornton (1977) for the shallow water condition: 


S(£,6) = Bg@£ -H(8)T(W_) (15) 


where: 


H(5) = (=) (cos (6)) (16) 
TT 
accounts for angular spreading, 


2w°r (Ww) = 
n nN 


} iy) 


rw) = © *(w){1 + — 3 —2 
sinh [2w_z(W) 


Ie, 


and 


g ily 
ie Ne (18) 
where 
(Ww) 
satisfies the condition 
z Tanh (WZ) = 1 (195 
(Ww) = 1 in deep water, and for W, less than l, 
. ~ ee? 
r (We) = x Wh (20) 
and 
il -3 
S(f,6) = zagdf H (96) (2 
where q is a coefficient on the order of Las a value of 


.073 is used in Wang's model. 

The bottom dissipation term, Ege is due to work done on 
the bottom by the bed shear stress. The bed shear stress is 
proportional to velocity squared and the rate of energy 
dissipation is proportional to the cube of the velocity. 


Hasselman and Collins (1968) proposed a quasi-linearized 
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expression for a frequency-de ident dissipation function: 
Sey 2) aa c,|U, | up (£) (22) 


where Ce is a bed shear stress coefficient, u, (£) is the wave 
induced bottom velocity at a particular spectral component, 


and U, is the bottom velocity for the total flow field. 


b 
Applying linear wave theory transfer functions relation the 
energy spectrum of the surface elevation to the velocity 


field, (22) can be expressed as: 


2,2 
a solace aiih> MRLs, cactle 
Eg (f) = pcel 1S (£)U, (23) 


a cosh (kd) 


where the term in the brackets is the spectral transfer 
function. The total bottom velocity is obtained by applying 
the spectral transfer function to the entire spectrum: 
MU) = (f =e se “k*s(e) af} (24) 
0 go cosh’ kd 
DD; THE SURF MODEL 

This model describes the surf zone width, number of lines 
of breakers, breaker type, breaker angle, significant and 
maximum breaker height, and the effective surf. The calcula- 
tions are based on empirical formula given in the Shore Pro- 
tection Manual (U.S. Corps of Engineers, 1977). The input 


wave 1S monochromatic, and requires only significant wave 


Pall 


height and period. Initial wave height and direction are 
determined for the entire model grid using swell direction 

and wave height calculations described in Chapter II, sections 
Band C. The depth of the first breaker is determined by 
comparing the calculated wave height with the incipient break- 


ing criteria (Noda et al., 1974; Thornton and Guza7 ei aa 


Ho 


ao 
= 0.43 pant) (25) 


The width of the surf zone, Xs is then calculated using 
the relationship between breaker height, breaker depth, wave 


period, and bottom slope (Weggel, 1972). 


Hy (x) Hy (x) 
a, Uy = b(n) > ae an (26) 
T 
assuming 
Hp 6) = ct dX, (26a) 
dy (x) = ert £X) (26b) 
Solve: 
(a'b£) x, + (ated + a'cf - b'f + d)X, + (a'ce + c-b'e) = 0 
(26c) 
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The coefficients c, d, e and f are determined by uSing 
standard slope intercept techniques. The positive root of 


this quadratic equation is the width of the surf zone. 


a(m) = 1.36(1-e 7°") (27) 
b(m) = aaa = (28) 
l1+e : 
a = alm) 
a = 
oe 
At 
b' = b(m) 
where: 
m = bottom slope between grid points 


The first breaking wave height, denoted as HBK in the 
program, is then determined using the surf zone width by 


interpolating between grids. 


JaGi=) inet Nee 


Te (Hee (2-0) (HT 


I,J slg 


eZ2) 


where H is the original wave height from energy considerations 
only and I and I-l are the grids that bracket the deepest 
breaker in the on-offshore direction. 

The depth of the first breaker, DBK, is calculated also 


using the surf zone width an c<he interpolation between grids: 
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e Re 


= Pa Ie 
DBK Wee eed gy) (4-0) + (Dry 3) (2 iL) ieee eG ) xX 


Jl 


(30) 


where D is the bottom depth in the program. 
Additional lines of breakers shoreward of the initial 


breaker are calculated by: 


t 
y = p dy, next (31) 

b next l1l+aitd 

b next 
where 
h -h 
i I+l,J Aree) 

qd a oa Neg fs —— io ae (I-1) Ax] (32) 


The angle of the breaking wave (Angz) is calculated by 
interpolating the distance between grid points and the linear 
change in the wave direction: 

poy eee 


Angz = Ang, . + ($19 49) (x - (1-1) Ax) (33) 


where Ang is the calculated wave direction at the grid point. 


The type of breaker is based on the surf parameter: 


= Mt (34) 


where: 
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m = bottom slope 
Ho = deep water wave height 
L, - isons deep water wave length 


Breaker type is specified by the following criteria: 


aa 2 surging or collapsing 
0.4 < ge, < Ze plunging 
ee 0.4 spilling 


Effective surf is calculated in accordance with 
COMNAVSURFPAC/COMNAVSURFLANT Instruction 3840.1 for use in 
naval amphibious operations. "Effective surf height is a 
numerical value, expressed in feet, and provides a guide for 
judging the feasibility of landing operations for any given 
landing craft under existing or forecast conditions." Effec- 
mave Surf is a function of: 

l. Breaker Height 

2. Breaker Period 

3. Breaker Type 

4. Littoral Current 

eee Wind direction and velocity 

6. Secondary wave height 

7. Breaker angle 
The parameters for the empirical calculation of effective 
Surf are calculated in the surf condition part of the model 
except for littoral current. The nearshore current circulation 


in the model is described in the next section. 
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E. NEARSHORE CIRCULATION MODEL 

The nearshore circulation is generated by changes in momen- 
tum due to breaking waves within the surf zone and by wind 
induced surface stress. The nearshore circulation model is 
based on the depth integrated and time averaged equations of 
motion by Phillips (1977). The breaking wave-induced currents 
are calculated by alternately solving the momentum and the 
continuity equations. The equations are solved numerically 
with the finite differences scheme taken directly from Kirby 
and Dalrymple (1982) which is an updated version of Birkermeier 
and Dalrymple (1977). The x and y momentum equations are 
integrated over the depth and time averaged over one or several 
wave periods. The equations contain bottom shear stress, mean 
surface shear stress and excess mean momentum stress due to 
wave action. To simplify the model, the formulation neglects 
the nonlinear mean momentum flux terms and the lateral (hori- 
zontal) transfer of turbulent momentum. The continuity equation 


1S given by: 


male 
ct|> | 


a ~ 3 = ~ i 
+ gylu(d + nd] + SF lvid +m] = 0 (35) 


where n is the free surface elevation, d is the still water 
depth, (u,v) are the (x,y) horizontal velocity components and 
the overbar indicates time averaging. The horizontal momentum 


equations are 





= = aS aS 
du .) 1 XX x = = 


0 (d+n) 


22 oe ve) gn (2 + SY - T+ T (37) 
0 (d+n) 
The oe AO; SES are the mean momentum flux due to the 
wave motion (commonly referred to as radiation stresses). 
es Tey are the surface shear stresses due to wind, and 
it y are the bottom stresses. 


The radiation stresses are given by (Longuet-Higgins and 


Stewart, 1964): 


Be = E[(2n - Wacos-6 + {(n - W/2) san 6 | (38) 
se Z 

Soy = E[(2n - 1/2)sin’8 + (n - 1/2)cos’ 8] (39) 

eee So = nein(ze) (40) 

<7 Ge Ae Sin 


where E iS wave energy = sOGH’ ; n= fil & 


The surface wind shear stresses are calculated (Van Dorn, 


1953) by: 


ae oK|W|W, (41) 


K|W|W 42 
ey oK|w| 6 | (42) 


where Woe Mh are the wind speed in the (x,y) directions, and 


K is the wind stress coefficient of order moe 


2 


The bottom shear stresses are linearized by assuming weak 


mean currents (U,V) (LaBlond and Mysak, 1978). 
4 
T = 9 -c, u_ U (43) 
T 
Se 3S (44 
prs ) 


where Cr is the Darcy-Weisbach bottom friction factor and 


Un, 1s the maximum wave orbital velocity at the bed: 


: Gla 
OO es sii ee (45) 


The maximum time step At is based on the linear stability 


criterion (Kirby and Dalrymple, 129382). 
/ 2 2 
Bee Fe (CAx)” + (hy) 97 Zag) (46) 


However, in practice, Kirby and Dalrymple (1982) recommend 
using a quarter of this value. The actual At used in the 


program is 


WZ 


2 2 
en = = a (47) 


(gd) 


The finite difference form of (36) and (37) are given 


by: 
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er > es 


( 4Ayod 4Ayod 
E + T E + T 
= SX1 5 SXq-1,J, DXr J P¥I-1,) ; 
2od 26d 
where: 
d +d 
qd = (eth is the average depth in 
the X-direction 
w = (My = We 
Bows SOS esi | Tey. 
S = S => 
eee es, Jo1) +1, J “*I-1,3) 
4Axod 4Axod 
= Tee whe 
ee Ease: : I,J) (@iaay york 
Ayod 2ed 
T + T 
Te ea Eel, 5 (49) 
2od 
where: 
d wc gs =2 
ah I,J-l1 
2 
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The finite difference form of (35) Solving formers 


— = wie I-l1,J-l 
= n + At(U (— : ) 
BOS oes jo 2AxX 
t+] 0 = a I,J° chy 
d ~~ 
1a) IE pal 
ss a ef ee 


The new water depth is then given by: 


The nearshore circulation solutions are "Spun up" to 
steady-state condition. The model is started from state of 
rest with zero wave field present (Kirby and Dalrymple, 1982). 
To reduce the effect of seiching or instability, the wave 
height H at the offshore grid is gradually brought up <o 


its full value uSing: 


7 Z0 
H = A tanh (=-) (545 
S 
where: 
t = model time 
aye = arbitrary fixed period 
Hy = full height of boundary wave 


30 


The hyperbolic tangent function is used to control the 
growt.. of the boundary wave during an iterative process. The 
quantity (2t/T_] increases linearly while the tanh function 


increases smoothly toward unity. 
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III. MODEL APPLICATION ON THE HEWLETT-PACKARD 9845B 


A. CONVERSION INTO BASIC 

The U.S. Navy contracted for the development of the wave 
and surf model which was written in FORTRAN. In its original 
form, the model was written for a large non-portable computer. 
To make the model useful tactically for mobile U.S. Navy forces, 
the program is rewritten specifically for deployed micro 
computers like the HP-9845B. It is required that the model 
be written ina user friendly format with common naval conven- 
tions and terms. Error checking and input review had to be 
programmed into the model to be operationally acceptable to 
the U.S. Navy. The Wave and Surf Model had to run on the 


HP-9845B to give reliable, useful and easily understood output. 


B. MODIFICATIONS IN THE WAVE AND SURF MODEL FOR THE HP-9845B 
The original contractor delivered a FORTRAN version of 
the wave and surf model written ina format of "stand alone 
Subroutines." This version was suitable for conversion to 
the HP-9845B because its size was within the micro computer's 
memory requirements. The subroutine format also allowed for 
easier testing of the model. The main change made in the 
program involved the use of common statements. The FORTRAN 
version made extensive use of labeled common blocks. The HP 


basic only allows one common statement. This minor problem 
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was solved by increasing the size of variables in the call 
statement to most subroutines. Only the most -rten used 
general variables were used in the single common statement. 
The language of the input statements is changed to agree 
with naval conventions. Wind and swell are described by the 
direction they arrive from. The beach is described from the 
point of view of a ship offshore. 
The interactive subroutines are expanded with lines of 
code that: 
1) explain the input needed 
2) ask for the data input 


3) error check the data input before accepting or 
rejecting and requesting correct input. 


The program was expanded in the calculation code with 
explicitly assigned variables for counters and matrices. For 
example, effective surf was stored into a unique array rather 
than sharing an array used for several other calculations. 
Empirical values from the COMNAVSURFPAC/COMNAVSURFLANT Instruc- 
tion 3840.1 were put into an array rather than use the FORTRAN 
version data fitting calculation. The additional memory 
space for these arrays was traded for explicitly assigned 
data and clearer programming. 

Labels were placed before every significant calculation 
in the model code. Spaces were added between code to separate 
functional groups. Indentations of the code were made to show 
loops and double loops. These changes added to the readability 


of the program and made the model logic easier to follow. 
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In addition, changes/corrections were made to the sea 
surf and swell model for undefined constants, incorrectly 
dimensioned arrays, logic errors and improperly defined 
angles. Some of these changes were made with the cooperation 


of the contractor. In the program code, each change is 


documented. 


a 


iv. a wtoCUSSION 


A. USABILITY 

The sea, surf and swell model is programed to be inter- 
active in language a sailor can understand. The units for 
wave height and bottom depth entires can be either meters or 
feet. Grid spacing can be in nautical miles only for the 
open water version. However, the nearshore model allows 
either feet or meters. All entries are requested with a 
format example. Where possible, the format example is in 
standard World Meteorological Organization code, which should 
be familiar to all U.S. Navy Aerographers Mates. Suggested 
input values, or default values, are offered where appropriate 


to make the model more "user friendly." 


B. DATA INPUT REQUIREMENTS 

The quality of the input data is the most important factor 
Peeemmining accurate wave and surf forecasts. The bathymetry 
begins to affect sea and swell when the depth is less than half 
a wave length. For the same depth of water, the effect of 
the bottom will be strongest on longer period waves. The 
more shallow the water, the greater the effect will be on all 
waves. Accurate bathymetric data is needed in deeper sections 
of the model grid, but becomes increasingly important for 


shallow waters. The wave period and depth is used to compute 


BES 


wavelength and phase speed. These parameters determine 
forecast wave direction, energy, surf type and currents 
generated in the surf zone. For the same reason, initial 
wave direction is as important as the wind direction. The 
bottom type input influences energy dissipation and surf zone 
currents. The surf model calculation is least sensitive to 
the bottom type. 

The most difficult, and clearly the most important inoue 
is the offshore wave input. Careful consideration must be 
given to the many sources of deep water waves. The open 
model sums the energy from all deep water wave inputs. The 
wave spectrum input must be accurately described. The ocean 
can be nearly flat concealing long period swell that can 
generate catastrophic surf on beaches. The ocean can have 
fully developed seas generated by continuous wind. Ina 
fully developed sea, the energy of the waves is distributed 
over a wide range of periods that must be entered into the 
model for accurate surf forecasts. Only with a sound knowledge 
of deep water wave conditions can the model be used to fore- 
cast shallow water conditions. 

The open water part of the model has the most options for 
wave energy input, which include: 

1. A single narrow band swell. 


2. A significant wave height and period which can be 
expanded in terms of a Bretschneider spectrum. 


3. Specific energy levels can be entered at appropriate 
frequencies such as from the Navy Spectral Ocean Wave 
Model. 
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It is important to note that the output from this model is 
a Single significant wave heicnt and period moving in the 
direction of maximum energy. This could be a poor represen- 
tation if the seas were large and from diverse directions. 
The surf and nearshore circulation parts of the model 
allow only a monochromatic wave input. This input could be 
directly from the open water calculation or an input of the 
user's choice. If there 1s more than one dominant wave group 
approaching the beach, the model representation of the surf 


zone will be incomplete. 


See OELECTING THE GRID SIZE 

The grid size and the distance between grid points will 
determine the model's ability to resolve the wave climate. 
The open water model has grid distances in nautical miles. 
The grid mesh should cover the area from deep water waves 
into shallow water, outside the surf zone. The grid size for 
the surf and nearshore is measured in feet, and the selection 
1s critical to resolving the surf zone. Kirby and Dalrymple 
(1982) outline the following requirements for the nearshore 
erreulation grid. 

1. The grid must extend offshore far enough to remove the 
offshore region of the domain from the influence of 
currents driven by the surf zone, and to allow for 
the specification of a uniform longshore depth which 
will not significantly alter the wave refraction 


results in the nearshore. 


2. The grid mesh must be fine enough to resolve the surf 
zone adequately. 


3. In the event that the effect of a single physical 
feature isolated in the longshore direction is to be 


a 


modeled, the longshore extent of the grid must be 

sufficiently large to isolate the physical system 

from the effect of images created by the longshore 

periodicity requirement. 

Experience running the model suggests that a minimum of 

3 grid points must be inside the surf zone. This would not 
be possible for waves breaking directly on the beach. A 
further constraint is the number of grids in this model (10 


offshore rows and 13 columns alongshore) due to the limita- 


tions of speed and storage area of the micro computer. 
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V. MODEL TEST AND EVALUATION 
A. DEEP WATER WIND GENERATED WAVES COMPARED WITH JONSWAP 

CURVES 

The model approximates the area under the spectral density 
curve (wave energy) with rectangles. The more rectangles 
used, the better the approximation will be, but computational 
time also increases. The open water model requires approxi- 
mately 45 minutes to calculate each frequency band on the 
HP-9845B. Narrow band swell waves might be closely approximated 
with a single energy band and a band width of .01 Hz. A 
broad band wave spectrum, such as wind waves, needs to be 
represented by multiple energy bands. 

To test the wind generation model and the effects of 
bottom friction and band width selection, the model was run 
for constant depths of 100, 10, and 5 meters (Figs. 3, 4, and 
Bee the initial wave height is zero. The bottom friction 
coefficient chosen is .0l for a fine sand bottom. The grid 
Spacing is 10 nautical miles resulting in a maximum fetch of 
90 nautical miles. The results are compared with the JONSWAP 
Significant wave height as a function of fetch from Bishop 
(1983). The model results favorably agree with the JONSWAP 
curve considering the relative coarseness of the energy 
bands. For the 100 m depths, bed friction was insignificant 


in changing the wave height as expected. 
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The model was run for the same input with the flat bottom 
feneimcwmanagea te LOM@meters™s The effect of friction only 
slightly reduced the fully arisen significant wave heights. 
However, the model generated fetch limited wave heights were 
approximately half a foot higher than the JONSWAP results. 

The increase in the height of the "shallow water waves" was 
probably a result of the change to a shallow water saturation 
spectrum. In shallow water the slope of the energy saturation 
Spectrum with respect to frequency is -3 (Thornton, 1977). 

In deep water, the slope is -5 (Phillips, 1966). This results 
in relatively more wave energy at higher frequencies for the 
shallow water spectra. 

To test the bottom friction in the model, a run was made 
with the same input as before except the depth was changed 
to 5 meters. The significant wave height curve shows a slow 
decrease after the first grid point. The model calculates 
larger frictional dissipation as the down grid distance in- 
creases. The model nearly reached a constant wave height, 
balancing generation and dissipation forces after 80 nautical 
miles. The values of predicted wave heights by the JONSWAP 
spectrum agree with the present model within 10 nautical miles, 
but tend to be larger for a longer fetch. For an open sea in 
a shallow region, the JONSWAP spectrum ignoring bottom friction 
effect gives a wave height four times larger than the pre- 
dicted by the present model; this result emphasizes the 
importance of including bottom friction in shallow water wave 


models. 
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The wind generation part of the model was again tested 
for 30 kt winds blowing over the same area. The significant 
wave heights computed by the model increase rapidly as the 
downwind fetch is longer. The model predicts a near fully 
arisen sea after only 50 nautical miles. The JONSWAP model 
predicts approximately linearly increasing wave height values 
for 30 kt winds until reaching fully arisen seas after 200 
nautical miles of fetch. The large difference between the 
model calculations at 30 kts and the JONSWAP curves suggest 
the model over-builds seas for this wind speed and grid spacing. 

The energy transport calculation in the model was found 
dependent on grid spacing and is invalid for large grid 
spaces and higher frequencies. Negative energies resulted for 
input at frequencies greater than .3 Hz with a grid spacing 
of 10 nautical miles. By decreasing Ax = Ay to one nautical 
mile, the model can compute energy transport for .4 Hz and 
higher. A "fix up" in the program is made by setting any 
negative energy values to zero. The inability of the model 
to compute energy transformation for the short waves should 
be noted, even though this is the low energy side of the 


spectra. 


B. NEARSHORE CIRCULATION ON A SLOPING BEACH 

The nearshore current prediction part of the model was 
tested on a uniform sloping beach. Since the depth is uni- 
form in the alongshore (Y) direction, the wave height, wave 


direction, current, wave set-up and surf should also be 
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Figure 6. 


uniform in the Y direction. The calculated wave directions 
were constant in the Y direction, the wave height and surf 
Changed a few percent but the variations of current and wave 
set-up in the Y direction deviate more than 20 percent. 
The wave height and surf were slightly inconsistent with ex- 
pected results and this is thought to be caused by the finite 
differencing scheme producing errors from the boundaries. The 
larger variation in the current and wave set-up values is 
probably due to the inconsistent mass flux calculations. The 
boundary errors induced in the wave calculations resulted in 
a nonuniform forcing term, and, therefore, they may not be 
able to describe a constant current. The model is expected to 
give best results in the interior of the computational domain 
for wave and surf prediction. 

The nearshore current calculation portion of the model 
was compared with data acquired at Santa Barbara on 4 February 
1980. Initially the currents calculated were one order of 
magnitude smaller than the observations. 

Consequently, three basic problems with the surf model 
were identified: 

1. The model computed only 36 seconds of surf after 
beginning with still water. The contractor's program 
calculated only 10 time steps, or 3.6 seconds. This 
length of time was too short to even approach a 
Steady state current. Kirby and Dalrymple (1982) 
always ran their nearly identical version of this 
model for a 1000 time steps representing more than 
5 minutes real time. 

2. The resolution was too coarse to describe the current 


field. The finite difference scheme is unable to 
compute currents nearer than 2 gird spaces from a 
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depth of zero (at the beach presumably). The surf 
zone for this test was closer to the beach than 2Ax 
or 20 meters. 

3. The changes in radiation stress, which forces the 
current, were very small across the breaker line where 
it should be the largest. The current would remain 
too small even with items 1 and 2 corrected due to 
weak forcing. Realistic radiation stress was not 
being calculated because the model gave incorrect 
wave heights for breaking in the surf zone. This 
correction was made. 

The model was run again for comparison with Santa Barbara 
data. The grid size was reduced to Ax = 2.5 meters and 
Ay = 5m The time step for this small grid was approximately 
-6 seconds. The HP-9845B was allowed to calculate for 381 
time steps which took 9 hours. This calculation represented 
5.28 minutes real time for waves entering still water. The 


maximum stable time step is a function of the grid size 


TMAX = 1/3 a 
The current was approaching a steady state with a very slow 
rate of convergence (Fig. 7). It would take more than 23 
hours of CPU time on the HP-9845B to run the model for a 
1000 steps as done by Kirby and Dalrymple (1983). Attempts 
to increase the time step and decrease the CPU time resulted 
in numerical instability (Figure 8). The model was run uSing 
the same input except the At was increased 5 percent each 
time step until reaching a maximum of one second. Instability 
began to occur after 20 seconds with values oscillating around 


the stable solutions. After 25 seconds the unstable model 
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Figure 8. 


The stable 


California. 
increase in current is shown with the 


off Santa Barbara, 


The unstable calculation is 


shown with the heavy dashed line. 


solid line. 
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of the current became large and changed sign. The surf model 
is very slow to converge to a solution and unstable for even 


slightly longer time steps. 


Cw SURF COhberron 

The surf zone wave heights and directions calculated by 
the model were compared with observations. The comparisons 
were made with observations for two California beaches at 
Santa Barbara and Torrey Pines. These beaches are nearly 
plain sloping. The offshore boundary condition of wave 
height, period and direction are given in Table I. The 
beach material for the model friction coefficient selection 
was fine sand (Cf = .01). The modeled and observed wave 
direction agreed, within round-off, to the nearest degree 
for all days. The waves turned toward normal while crossing 
the surf zone exactly as expected. 

Model comparisons for observations on 3 February are shown 
in Figure (9) which shows the model slightly under-predicts 
the wave heights but follows the general change across the 
surf zone. The breaker line was forecast to be 14.9 meters 
from the shore with a height of .68 meters. The observed 
breakers began more than 40 meters from the beach with a 
height of .68 meters. The model incorrectly forecast one 
line of breakers when at least two were observed. 

The 4 February results (Figure (10)) show reasonable 
prediction of wave height. But using the wave height 


information, the model did not predict the location or height 
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TABLE I 


Offshore Boundary Conditions 


Date Depth es - Wave direction 


(m) (m) (sec) ig) 


Santa Barbara, Leadbetter Beach, California 


3 Feb 80 4.1 Zo 14.3 =f sis: 
4 Feb 80 4.2 56 14.3 =3720 
5 Feb 80 4.1 45 TZ -8.4 


Torrey Pines Beach, California 


4 Nov 78 6.7 235 14.3 0 


10 Nov 78 6.6 .68 152.9 0 


51 


30.0 





SERRE MnEtaRS 


OPSTPNGE PRM ol rie 


Model computed values are 


represented by the solid line, observed 


Santa Barbara surf zone wave height on 3 
heights are dashed. 


February 1980. 


Figumes 9. 


oz 


ee ee i ee a OO OO | 


a -— =  . =< : 
’ ' . ‘ 

’ | ' 

t ‘ t 

4 ‘ ' 


| 4 1 

eg eee oS an ae 
'‘ ‘ 

| 

\ 


$0.9 


me a 


mse = =e = Oe eae ne tee a 


1 ‘ 
f ' 
1 i 
1 i] 
‘ i] 
1 t] 
' 4 
i] 

1 4 
= Kes u 
‘ 

4 

1 

i] 

t 

' 

4 

! 

! 

4 

4 

' 

J 

t I 
1 ' 
1 t 
' ‘ 
t] ‘ 
! 1 
‘ t] 
' ' 
' ' 
Beret k 1 
| 

4 

' 

t} 

' 

t 

' 

L] 

60.9 7f 
FFOIM TRE BEACH [ts NETTES 


| 
| 
| 
| 


| 






ah a) 
(eee eee 


pe eee Gee 


20.0 











t 

1 ; { oO 
| Se iat ata eo Bee _. - epee ee ere a a Se +S 
Bea ; | ! 
| | 4 4 1 
: : 
’ ' = 

—— ae — 2 ee eS ! 
O°} 8° 3°0 vaG cae ‘omre) 

| GUSLGW NI LEST GARM | 
| | 
| 
| Se ee eee ee Se 


Figure 10. 


ght on 4 


Model computed values are 


Santa Barbara surf zone wave hei 


February 1980. 


observed 


? 


represented by the solid line 


heights are dashed. 


D5 


or breaker. The observation indicated two lines of breakers, 
the first near 60 meters and the second near 38 meters. The 
maximum measured wave height was .67 meters occurring 60 
meters from the beach. The model forecast only one line of 
breakers .84 meters high and 16 meters from the beach. 

The 5 February results (Fig. (11l)) show the model almost 
forecast the surf height exactly at grid points. The model 
forecast one line of breakers 11 meters from the shore with 
a height of .6l1 meters. The observed breakers began about 
39 meters from shore with a height of nearly .60 meters. 

The model was also compared with surf zone observations 
from Torrey Pines State Beach, California for 4 and 10 
November 1978 (Table I). The average angle was near normal 
to the shore. Fig. (12) shows the model input boundary wave 
on 4 November height was too small by about .05 meters (2 
inches). This small error was due to interpolation between 
sensor locations. The model followed closely the observed 
wave heights between 160 and 240 meters, but it over-built 
the wave height inside 160 meters. The model placed the line 
of breakers 82 meters from shore with a height of .76 meters. 
Observed breakers were near 160 meters with a height of .48 
NeEEeES. 

The 10 November 1978 model of Torrey Pines was based on 
boundary conditions entered at 6.6 meters. Figure (13) shows 
uniform over-prediction of the wave height at grid points. 


The model forecast two lines of breakers, the first line 133 
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meters from the beach with a height of 1.7 meters and the 
second 83 meters from the beach with a height of .75 meters. 
However, the obServations showed no corresponding lines of 
breakers. The observed heights of the lines of breakers are 
./2 meters or less. 

The wave height calculation by the model is, in general, 
close to the observed wave height. The surf zone width 
based on (26) iS incorrect in every test run. The resulting 
calculations of breaker position, angle and height are 
dependent on the width of the surf zone and therefore these 


forecasts are incorrect. 
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VI. SUMMARY AND CONCLUSIONS 


A. SUMMARY 

A numerical model describing surface gravity wave trans- 
formation in shallow water region is written in BASIC for use 
on the HP-9845B. The original version of this model was 
based on the work by Wang and Chen (1983). The model consists 
of an open water wave program and a surf program. The open 
water model considers a wave energy spectrum generated by 
winds; the surf model predicts a shallow water wave field and 
the shear currents induced by wave breaking on an arbitrary 
bottom topography. The model is tested for cases of constant 
depth, a uniformly sloping beach and actual beaches. The 
application of the model to a field beach indicates the poten- 
tial value for naval operational planning. Combinations of 
different conditions of wind, wave, and bottom bathymetry can 
be effectively simulated and serve as the basis for determin- 
ing an amphibious landing. 

The open water model is used to forecast sea heights over 
continental shelf and wide stretches of shallow water coast. 
The calculations of wave energy are the result of wind forcing, 
Peet eaistance, and bottom friction. The model predicts the 
spectral distribution of wave energy. 

The surf model can predict wave shoaling, coastal cur- 


rents and multiple breaker lines. In the model, the wave 


5) 


heights and directions are calculated for a section of beach. 
The breaker heights, breaker types, and location of severe 
breakers are empirically determined based on the wave height 
calculations. This provides a map of the surf zone in two 
dimensions for Navy operations. 

One of the drawbacks of this model is that the model is 
not programmed in an efficient manner. The wavenumber at a 
local water depth is calculated over a hundred times in nearly 
every subroutine for different frequency components. To allow 
for various options chosen and setting up the connectivity 
between subprograms, the integrated model was too large and 
not very well structured. Moreover, the unsteady approach 
to the steady state solution is time consuming and needs 
Special care to avoid numerical instability. Using a steady 
state approach and reforming the program structure would 
facilitate the application of this model and improve the 
operational efficiency. Further modeling efforts should be 
directed toward this direction. 

A problem in running the model is the choice of the 
energy bandwidth. The open water model uses energy spectra 
to describe the wave field. The selection of energy bands 
and bandwidths requires knowledge of wave spectra and gravity 
wave theory. The alternative to user choice is to have the 
model calculate energy contribution over the entire wave 
Spectra. This would ensure all energy contributions across 
the spectra are considered. However, CPU time would increase 
to more than 6 hours if only 10 energy bands were used. (The 


Navy Spectral Ocean Wave Model uses 15 energy bands.) 
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B. CONCLUSIONS 

This sea, surf and swell model at present is not suitable 
for operational use. The separate regimes, open water and 
nearshore current and surf, all have problems. 

The open water part of the model requires sophisticated 
spectral energy input. The average Aerographer Mate does not 
have the knowledge to properly select input parameters. Only 
when the sea is dominated by a single narrow band swell with 
light wind is the wave energy input fairly easy (input wave 
height, direction and period). In addition, the open water 
model was shown to build up the sea height too fast for wind 
speeds of 30 kts raising doubt concerning the wind energy 
transfer part of the model. 

The nearshore current part of the model is very slow to 
converge requiring more than 23 hours to complete the neces- 
sary 1000 iterations. The numerical scheme is therefore 
unsatisfactory for the HP-9845B computer. Recommendations 
include: 1) using a more stable scheme to increase the time 
step and reduce the number of iterations; 2) incorporate an 
implicit numerical scheme to calculate a direct solution of 
the currents such as in Wu and Liu (1985); or 3) use a simple 
one-dimensional current model implying a uniform beach 
alongshore. 

The surf part of the model is probably the most "user 
friendly" and has the most potential. The surf part of the 


model is fast running and gives good solutions for grid point 
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wave heights. However, the width of the surf zone calculation 
1s very poorly predicted. This error in turn causes the 
breaker line position, height and angle to be incorrect. 

The incorrect breaker predictions cause the effective surf 


calculation also to be invalid. 
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